Matter wave quantum dots (anti-dots) in ultracold atomic Bose-Fermi mixtures 
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The properties of ultracold atomic Bose-Fermi mixtures in external potentials are investigated 
and the existence of gap solitons of Bose-Fermi mixtures in optical lattices demonstrated. Using 
a self-consistent approach we compute the energy spectrum and show that gap solitons can be 
viewed as matter wave realizations of quantum dots (anti-dots) with the bosonic density playing 
the role of trapping (expulsive) potential for the fermions. The fermionic states trapped in the 
condensate are shown to be at the bottom of the Fermi sea and therefore well protected from thermal 
decoherence. Energy levels, filling factors and parameters dependence of gap soliton quantum dots 
are also calculated both numerically and analytically. 



INTRODUCTION 

The possibility to cool atomic Fermi gases below the 
Fermi temperature Q with the developement of the sym- 
pathetic cooling technique has raised a significant 
amount of interest in theoretical and experimental stud- 
ies of atomic Fermi gases embedded in Bose- Einstein con- 
densates (BEC). Such systems have very interesting fea- 
tures which allow to explore a whole class of phenom- 
ena, ranging from strongly correlated effects in fermionic 
atomic gases to novel superfluid phases in Bose-Fermi 
mixtures (BFM) analogous to BCS superconductivity. 
The possibility to trap fermionic states into localized 
bosonic components of BFM makes these systems also 
of potential interest for quantum computing. The basic 
element of a quantum computer is the so called quantum 
bit, a special quantum state usually realized with a two- 
level system, which is manipulated by logical quantum 
gates (unitary operators) and which must be preserved 
intact for an extended period of time to allow compu- 
tation 3]. An important role in this context is played 
by the interaction of the quantum system with the envi- 
ronment (thermal reservoir) which unavoidably produces 
decoherence. For an ideal q-bit one would like to keep 
the quantum system well protected from decoherence but 
not completely isolated from the environment to allow 
gate manipulations. This condition could be realized in 
Bose-Fermi mixtures by trapping fermionic states into 
localized BEC, so to form quantum dots with two (or 
more) fermionic levels inside, and by using Feshbach res- 
onances to control the boson-fermion interaction, i.e. the 
coupling of the quantum bit with the environment. 

The aim of this paper is twofold. From one side, we 
investigate the possibility of using ultracold Bose-Fermi 
mixtures to form matter wave quantum dots with two 
or more fermionic levels inside. To this regard we con- 
sider solitonic BEC states as trapping (expulsive) dot 
(antidot) potentials and show that the fermionic bound 
states formed inside are at the bottom of the Fermi sea 
and therefore well protected by thermal decoherence. 

From the other side, we investigate the spectral prop- 



erties of a BFM both in parabolic and in periodic poten- 
tials and demonstrate the existence of BFM gap-solitons 
in optical lattices (OL). 

To this regard, we derive equations of motion for the 
mixture adopting a mean field approximation for the 
bosons and assuming non interacting (spin polarized) 
fermions in true quantum regime. These equations are 
solved by means of a self-consistent procedure and the 
spectral properties of BFM, both in parabolic traps and 
in optical lattices, derived. In particular, we show that 
fermionic bound state levels at the bottom of the Fermi 
sea can be trapped inside gap-solitons for both repul- 
sive and attractive boson-boson interactions. We also 
show that, as the boson-fermion interaction is varied, the 
trapped levels undergo avoid crossings or level crossings 
according to the symmetry of the corresponding wave- 
functions with respect to the OL. Analytical expressions 
for energy levels and for the number of fermionic bound 
states trapped inside attractive gap solitons are also de- 
rived and the results compared with those obtained from 
direct numerical self-consistent calculations. We find 
that the filling factor of a gap soliton follows the law 
\/\Xbf\Nb, where iV^ is the number of bosons, nj 
the number of fermions trapped in the condensate and 
Xbf the strength of the boson-fermion interaction. The 
possibility to change Xbf in a wide interval of values using 
Feshbach resonances suggests the possibility of an exper- 
imental check of this results. In particular we expect 
that, as the boson-fermion interaction is increased, the 
progressive filling of the gap soliton should give rise to a 
depletion of fermions in the OL which follows the above 
law and which could be monitored by imaging techniques. 

The paper is organized as follow. In section II we in- 
troduce the model and derive quasi-classical equation of 
motion of Bose-Fermi mixtures directly from a second 
quantized formulation of problem. In section III we in- 
vestigate the energy spectrum for BFM in attractive and 
repulsive parabolic potentials and discuss quantum dots 
(antidots) analogies. In section IV we consider BFM in 
optical lattices and show the existence of gap solitons 
both for repulsive and attractive boson-boson interac- 
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tions. The number of fermions trapped in a BFM gap 
sohton and the energies of the bound states are also cal- 
culated both numerically and analytically. Finally, in 
the last section the main results of the paper are briefly 
discussed and summarized. 



EQUATIONS OF MOTION FOR A BOSE-FERMI 
MIXTURE 

To derive equations of motion of a BFM it is convenient 
to start from the second quantized Hamiltonian of a set 
of bosons and fermions in interaction 
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are single particle hamiltonians accounting for the ki- 
netic energy and for the external potential confining the 
bosons and the fermions, respectively, while ?/j^(r), '!/)f,(r), 
'(/'j-(r), 7/;/(r) are usual bosonic and fermionic field op- 
erators satisfying [tpb{r),'4'b{r')] = [ipl{r),'4'l{r')] = 
0, [V'b(?'), = H''" ~ ''"') (and similar relations for 

the fermionic fields with commutator replaced by the 
anti-commutator). Notice that the fermions are assumed 
spin polarized and therefore non interacting, so that only 
the boson-boson and boson-fermion interactions appear 
in the Hamiltonian. Using the above operator algebra, 
the Heisenberg equation of motion for the bosonic and 
fermionic fields can be written as 
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ih^ = + j dr'i^l{r')Ubf{r' - r)i^f{r')4>f{r). 

For dilute mixtures one can replace the interaction po- 
tentials with the effective interactions 
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r-r'), Ubf{r-r') = gbfS{r-r'), (2) 



where the coupling constants gb,gbf, are related to the 
s-wave scattering lengths ab,abf, by the relations: gb — 
A-Kfi^ab/mb, gbf — 2'Kh?abf /mbf, where nibf denotes the 
reduced mass nibf — mbmf/{mb + ruf). Using these 
potentials, the above equations reduce to 
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where pb = ^lipb and pf = tp^ftpf, are the bosonic 
and fermionic particle densities, respectively. Expand- 
ing Bose and Fermi fields in terms of the usual creation 
and annihilation operators a^, d and c^, c, for bosons and 
fermions, respectively, 

and replacing all field operators appearing in the bosonic 
equation with their expectation values on the ground 
state (mean field approximation) 

^fc ^ (A) ~ 4ir,t), Pb ^ (pb) « Nb\i;'oir,t)\^ 

Nf 

Pf^{Pf)^J2\^^r.ir,t)\', 



we have that Eqs. Q reduce to 
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The completeness of the fermionic Fock space implies 
that the second equation in Q can be satisfied only if 
the terms in the square bracket vanishes for all n , this 
giving the following Nf + 1 coupled equations 
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with Pf — IV'i^ P the fermionic density and ■00 = ■ 

From the above derivation it is clear that while bosons 
are treated in a classical mean field approximation the 
spin polarized fermions are in true quantum regime due 
to the lacking of interaction among them. Also note 
that for gbf = Q the above equations separate into the 
Gross-Pitaevskii equation for the BEG part and Nf iden- 
tical (due to the lack of interaction among fermions) 
Schrodinger equations for the fermionic part of the mix- 
ture. Equations (|3J) were also derived in 0, Q by means 
of a Lagrangian density approach. 

In the following we shall consider Eq. in a quasi 
one dimensional context by assuming the bosons and 
fermions confined in a cylindrical trap with transverse 
trapping frequency uj_\_ and negligible x-axial confine- 
ment. We also assume that the ground state energy of 
the transverse trapping potential hLO^_ is larger than the 
3D ground state energies of the bosons and the fermions. 
Within this approximation the bosonic and fermionic 
wavefunctions can be factorized as V = V'(^)V'-l(2/; -2^) 
with the transverse wavefunction taken as the ground 
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state of the trap potential independently of the longi- 
tudinal behavior and statistics [g. By integrating over 
transverse coordinates one obtains effective ID equations 
which are formally identical to Eqs. ^ but with boson- 
boson and boson-fermion interactions rescaled according 
to: gb ~> 2fuj_igb, gbf — > 2huj±gbf. 

Notice that the fact that fermionic and bosonic den- 
sities appear in the equations of bosons and fermions, 
respectively, suggests to consider them as effective po- 
tentials and solve the quasi-classical equation of motion 
by means of a self consistent procedure. Starting with 
a trial function for the bosonic part, we solve the quan- 
tum eigenvalue problems for the fermionic part, compute 
the fermionic density, derive the bosonic wavefunction by 
solving the bosonic equation and iterate the procedure 
until convergence is reached (an efficient implementation 
of this scheme for single BEC was discussed in 0)- In 
the following we apply this method for BFM in parabolic 
traps and in optical lattices, concentrating, for simplicity, 
on the quasi one dimensional case. 

BOSE-FERMI MIXTURES IN A PARABOLIC 
TRAP 

In this section we consider BFM in a trap poten- 
tial of the form Vtrap — ^muj'^{x — xq)^ and assume 
for simplicity rub = ruf, V^^^ = V/^f = Vtrap (the 
extension to the generic case is straightforward). We 
normalize Eqs. © by measuring space in harmonic 
oscillator length units ao — y^h/muj, time in units 
of h/Eo, with Eq the zero point energy of the oscil- 
lator. Introducing the parameter Xb(bf) gb(bf) — 
iSTTh^ujj_abg(bfo)/'m)Xb(bf), where abo,abf^ are scatter- 
ing lengths in zero magnetic field (we use Xb(bf) to 
change nonlinearities by means of external magnetic 
fields via Feshbach resonances), and rescaling wavefunc- 
tions according to ■0^ {4:TrgboNb/muja'^y/^^'' , 
{4:ngbfg/muja'^)^l}^, with a±_ = {h/muj±^)^^'^ , we obtain 
the normalized equations 

lijt = -vL + v^" + xfclV'T^' + XbfPfi^', 

= -^{.. + + «X&/IV'T] (6) 

where V = e{x — xq)'^ and a = abf/ab assumed in 
the following of the order of unity (notice that with 
this normalization the physical number of bosons Np is 
obtained from the dimcnsionlcss number Nb as Np ~ 
muia^Nb/^Trgbo). 

In Fig^ we depict the bosonic and fermionic densi- 
ties obtained with the self-consistent approach for a BFM 
with repulsive bosonic interactions {xb > 0) and attrac- 
tive bosons fermion interaction {xbf < 0). We see that 
while the bosons remain confined in the middle of the 
trap, the fermionic density is quite extended in space 
and has a hump corresponding to the bosonic compo- 
nent. This hump originates from the contribution of the 




FIG. 1: (a). Fermionic and bosonic densities for the case 
of Xb = IjXb/ = ~10, Nf — 20 and normalized number of 
bosons A^i, = 1, for a parabolic trap of strength e = 0.0025 
centered at xo = L/2 = 207r. (b) Effective potentials and 
energy levels of the Bose (open circles and dotted line) Fermi 
(filled circles and continuous line) mixture with densities given 
in panel(a) (horizontal line refers to the Fermi level). Panels 
(c) and (d) are the same as panels (a),(b), respectively, except 
for Xbf = 10. Plotted quantities are in normalized units. 



fermions trapped by the BEC acting on the fermions as a 
potential well (trapping impurity). Similar results were 
also reported in p. 

An interesting problem which has not been investi- 
gated is the distribution of the bosonic and fermionic 
levels in their effective potentials. This is done in Fig. 

from which we see the presence of two fermionic lev- 
els formed below the parabola at the boottom of the 
Fermi sea. Notice that the effective potential seen by the 
fermions deviates from the original parabolic one only in 
the center of the trap due the presence of the potential 
well created by the bosonic density. Thus, the effective 
potential seen by the bosons is much deeper and wider, 
due to the contribution of the wide fermionic density and 
of their self-interaction. From this figure it is clear that 
two fermionic states have lowered their energies by ad- 
justing themselves into the potential well created by the 
BEC giving rise in this way to a matter wave quantum 
dot i.e. a localized BEC with one or more fermionic levels 
trapped inside. Notice that the two level system exists at 
the bottom of the Fermi sea and are separated by a gap 
from other fermionic levels which rapidly became equally 
spaced, as one would expect for a parabolic trap. Also 
notice that the ground state energy of the bosonic part 
of the mixture is lowered inside the potential well corre- 
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FIG. 2: Time evolution of bosonic and fermionic densities in 
Fig. ^ (left panel) and Fig. 0; (right panel) obtained from 
the quasi-classical Eqs. J^l on a line of length L = AQ-k. Space 
and time are in normalized units. 



spending to the hump in the fermionic density and that 
excited states, although lowering their energies, remain 
almost equally spaced (this is due to the fact that the ef- 
fective potential created by the fermions, except for the 
hump in the middle, looks still parabolic). 

It is remarkable that this situation remains stable un- 
der time evolution, as one can see from the left panel of 
Fig. [3 In Fig. we show the situation for repulsive 
Bose-Fermi interactions (x&/ > 0). From Fig. ^ we see 
that in this case fermions are expelled from (instead of 
attracted to) the region where bosons are spatially con- 
centrated, this increasing their kinetic energy and ris- 
ing a pressure on the condensate. Also notice from Fig. 
^ that the corresponding energy levels are shifted up 
in energy, much more for bosons (which have an effec- 
tive potential very deformed at the bottom) than for the 
fermions (which see almost the original parabolic trap 
but with a narrow hump in the middle) . From this figure 
it is clear that the case of repulsive interactions (xb / > 0) 
is also of interest because, although not producing trap- 
ping states, it induces correlations among the fermions 
as is evident from the deviations from equally spaced 
levels observed in the region across the Fermi surface in 
Fig. QJi. This case can indeed be viewed as a matter 
wave anti-dot and could be used to investigate effects of 
strong correlations among fermions induced by the BEC 
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FIG. 3: Panel (a). On-site symmetric densities for the case of 
Xb = 1, Xi)/ = ~li ~ 20 and normalized number of bosons 
A^i, = 1 for an optical lattice of strength e = 3 (thick line refers 
to bosons). To identify the symmetry of the state the OL is 
reported with an amplitude scaled by a factor of 60. Panel 
(b). Effective potentials (only lower part) and energy levels 
of the Bose (thin line and open circles) and Fermi (thick line 
and filled circles) mixture corresponding to the densities in 
panel(a). Lowest bosonic band (shadow region with negative 
slope) and lowest fermionic band (shadow region with positive 
slope) have also been shown. The fermionic band is half filled. 
Panels (c) and (d): the same as panels (a),(b), but for the 
inter-site symmetric case with Nb — 1.47 . Plotted quantities 
are in normalized units. 



V = £C0s(2a::), where e ~ Vo/Er- With this normaliza- 
tion the relation between the dimensionless number of 
bosons Nb and the physical one is Np = {kLa\/ab)Nb 
(typical values are in the range of 10"^ — 10^ atoms, for 
a± « lO-^m, fci = 10'^m-^,ab = 1 50 • 10-^°m). In the 
following we show the existence of gap-solitons in BFM 
and suggest them as matter wave quantum dots (anti- 
dots) (the anti-dot case will be discussed elsewhere 8]). 



BOSE-FERMI MIXTURES IN OPTICAL 
LATTICES 

In this section we consider a trap potential of the form 
Vtrap = Vbcos(2fcix), as a model for an optical lattice. 
We assume mb = rnf = to, V^^^ — v£.^ — Vtrap and 
normalize Eqs. (jsj by measuring space in units of k^^, 
time in units of h/Er {Er = H^k]^/2m is the lattice recoil 
energy) and rescaling wavefunctions according to Tpb 
{gboNb/27rEras_Y^^,l;b, V'/ ^ {gbfj2nEraiy/^y/^)^f. 
This leads to the same normalized equations (jSJ but with 



To this regard we consider BFM with repulsive 
bosonic-bosonic Xb > and attractive bosonic-fermionic 
Xbf < interactions. In Fig. ^ the bosonic and 
fermionic densities corresponding to an onsite symmetric 
(OS) gap soliton, are depicted (notice that the soliton is 
symmetric around a minimum of the potential). We see 
that while the bosonic density is localized in almost a sin- 
gle potential well, the fermionic density is very extended 
in space and presents a hump in correspondence of the 
BEC. In Fig. the effective potentials, the sketch of the 
band structure and the energy levels of localized states 
for both bosons and fermions are shown. Notice that the 
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FIG. 4: Time evolution of the OS (left panel) and IS (right 
panel) BFM densities in Fig. |^ and Fig. obtained from 
direct numerical integrations of Eqs. on a line of length 
L = 40n. Space and time are in normalized units. 



energy of the OS gap soliton lies in the gap between the 
first two bands as expected for repulsive bosonic interac- 
tions, while the fermionic level responsible for the density 
hump in Fig. I^l, lies below the lowest energy band at the 
bottom of the Fermi sea. As for the parabolic case, gap 
soliton states with trapped fermionic levels inside can be 
seen a matter wave realization of quantum dots. 

More complicated gap soliton states (or arrays of gap 
solitons) can also be formed. This is shown in Fig. O; 
for the case of a two hump BEG condensate obtained 
as inter-site symmetric (IS) gap-soliton (notice that the 
state is symmetric around a maxima instead than a min- 
ima of the periodic potential). We see that also in this 
case the fermionic density increases in correspondence 
with the two bosonic humps and two fermionic levels are 
formed in the corresponding effective potential (see Fig. 
131). We have checked that both the OS and the IS modes 
are very stable under time evolution (see Fig. 0)). 

Gap- solitons with different symmetries with respect 
to the OL, such onsite-asymmetric (OA) and intersite- 
asymmetric (lA) , can also exist in BFM but they appear 
to be metastable under time evolution. In Fig. |2K we 
show a gap-soliton of type lA together with two fermionic 
states trapped in the condensate. In panel (b) of this 
figure we depict the time evolution of the bosonic and 
fermionic densities of the lA gap-soliton as obtained from 
direct numerical integration of Eqs. ©. We see that 
after a time t « 60 the lA state decays into the stable 
OS wavefunction depicted in Fig. Notice that the 
lA-OS transition also induces changes in the fermionic 
states. In particular we see that the first excited state 
asymptotically decays into the Bloch state at the bottom 
of the fermionic band (see top curve of FigEt), while the 
lowest energy fermionic state remains well localized inside 
the gap soliton (see the time evolution reported in panel 
(c)). The lA-OS transition in Fig. occurs therefore 
with the expulsion of one fermion from the condensate 
and with a change of symmetry of the wavefunction of 
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FIG. 5: (a) Wavefunction of the intersite-asymmetric (lA) 
gap-soliton (thick line) plotted together with the wavefunc- 
tions of the lowest and first excited fermionic states trapped 
in the condensate (second and first curve from the top, re- 
spectively). The periodic potential, scaled by a factor 10 for 
graphical convenience, is also shown with a thin line at the 
bottom. Parameter values are the same as in Fig. |3 with 
Nb — 1. (b) Time evolution of the bosonic and fermionic 
densities of the metastable lA gap-soliton obtained from the 
quasi-classical Eqs. on a line of length L = 407r. Panels 
(c). Same as panel (a) but for the OS state obtained after 
the transition. Panel (d). Time evolution of the wavefunc- 
tion of the lowest energy state inside the condensate. Plotted 
quantities are in normalized units. 



the other fermion from lA to OS type (see second curves 
from the top in panels (a) and (c)). Similar phenomena 
occur for other metastable gap-solitons states existing in 
higher energy gaps and for different sign combinations of 
the interatomic interactions. 

Also notice from Fig. [SJd that the lA-OS transition oc- 
curs with emission of matter from the condensate mean- 
ing that the final OS gap soliton state contains less mat- 
ter than the initial state. This implies that the corre- 
sponding effective potential for the fermions is weaker af- 
ter the transition and therefore less effective to hold a sec- 
ond bound state in the condensate. As the boson-fermion 
(attractive) interaction and the number of bosons in the 
gap soliton are increased, we find that more fermionic 
levels are able to enter the condensate. 

An interesting problem to investigate is how a mat- 
ter wave q-dot (gap soliton) gets progressively filled with 
fermions as the product XbfNt is increased. To this re- 
gard we consider OS bosonic gap solitons with attractive 
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boson-fermion interaction and with both repulsive and 
attractive boson-boson interactions. In the left panel 



energy level of the quasi degenerate pair becomes non 
degenerate after the interaction. 




FIG. 6: Left panel. Energy levels of the fermionic bound 
states inside the condensate as a function of Xbf, for pa- 
rameter values Xb ~ li^ = 'i,Nb ~ l,Nf = 20. The top 
horizontal line represents the bottom of the fermionic band 
while the insets are enlargement of rectangles a,b, where level 
crossing and avoid crossing occur, respectively. Right Panel. 
Wavefunctions of the gap soliton (thick line) and of trapped 
fermionic states (ordered by decreasing energy from the top) 
a-t Xbf = —30. Other parameters are fixed as in the right 
panel. The vertical dotted line and the periodic potential at 
the bottom (scaled by a factor 60) were drawn to identify 
the symmetry of the wavefunctions. A vertical offset among 
curves was introduced to avoid overlapping. All plotted quan- 
tities are in normalized units. 

of Fig. we show the energy of the fermionic bound 
states inside the condensate as a function of Xb f for the 
case of a repulsive BEC (xt > 0) with a fixed number of 
atoms. We see that the presence of an arbitrary small at- 
traction between bosons and fermions is enough to trap 
at least one fermion in the gap soliton (notice that the 
energy of the fermion is below the edge of the lowest 
fermionic band, shown as a horizontal line in the figure) . 
This is a consequence of the one dimensional nature of 
the effective potential (for BEC in higher dimensions a 
threshold in Xbf is expected for the entering of the first 
fermion). From the inset (a) of the figure we see that 
two fermionic levels of almost equal energy enter the gap 
soliton at Xbf ~ —3.5. Notice that these levels come from 
extended (Bloch) states inside the fermionic band and are 
pulled below the lower band edge by the attractive boson- 
fermion interaction. The inset also shows that the two 
levels undergo a level crossing after which they proceed 
together as quasi degenerated states until a collision with 
a fourth level occurs. Notice that the last level detaches 
from the bottom of the fermionic band at Xfc/ ~ —28.6 
and, in contrast with the previous case, it undergoes an 
avoid level crossing with the quasi degenerated states at 
Xbf ^ —30 (see inset (b)). Also notice that the avoid 
crossing occurs with an exchange of degeneracy, i.e. the 
newly entered level, originally non degenerate, becomes 
quasi degenerate after the interaction while the lowest 
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FIG. 7: Bosonic density for attractive boson-boson interaction 
in an OL of strength e = 0.5. The continuous curve refers to 
the potential in Eq. Qwith a — 0.66 while the dots represent 
self-consistent numerical results. Other parameters are fixed 
as Xb = -l,Xi>/ = -10, iVt = 4.9, Nf = 20. Right panel. 
Gap soliton (thick line) and trapped fermionic wavefunctions 
ordered by decreasing energy from the top. Parameters are 
fixed as in the left panel. The vertical dotted line and the 
periodic potential at the bottom (scaled by a factor of 4) are 
drawn to identify the symmetry of the wavefunctions. A ver- 
tical off'set among curves was introduced to avoid overlapping. 
All plotted quantities are in normalized units. 

This behavior is in agreement with a general theorem 
of quantum mechanics, holding for hamiltonians depend- 
ing on a parameter ( with eigenvalues functions of the 
parameter) according to which the level crossing occurs 
between levels of different symmetry and avoid crossing 
between states with the same symmetry. In the present 
case the wavefunctions of the fermionic levels inside the 
BEC can be either symmetric or asymmetric with respect 
to the OL site around which the gap soliton is centered. 
Moreover, the symmetry alternates between consecutive 
levels, the lowest state being always of OS type. This 
is seen from the right panel of Fig. |2| where we have 
depicted the gap sohton found at Xfe/ = —30 with the 
four fermionic states trapped inside. Notice that while 
the second and third bound states have different sym- 
metries, and therefore their levels can cross, the second 
and fourth excited wavefunctions have the same symme- 
try and therefore their levels must give rise to an avoid 
crossing when approaching each other. This general rule 
for the crossing of the fermionic levels inside a gap soli- 
ton is in full agreement with the self-consistent numerical 
calculations shown in insets (a)-(b) of Fig. 

Analytical estimates for the energy levels and for the 
number of fermions in the condensate can be obtained 
for the case of an attractive BEC in shallow OL. To this 
regard we remark that for repulsive boson-boson inter- 
actions the bosonic density cannot be localized into a 
single potential well (due to the tunneling of matter into 
adjacent wells induced by the repulsive interaction), so 
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that satellites around a main peak always appear (see 
the right panel of Fig. 0) . In this case the corresponding 
effective potential in the fermionic Schrodinger equation 
gives rise to a spectral problem which is analytically dif- 
ficult to solve. For attractive boson-boson interactions, 
however, the bosonic density may have no satellites peaks 
(due to the attraction) and can be easily approximated 
by integrable potentials (this is especially true when the 
OL is shallow). 

In the following we consider the case of an attractive 
BEC with Xb = —1 in an OL of strength e = 0.5 and 
approximate the bosonic density with a Poschl- Teller po- 
tential (soliton potential) of the form 



Pb 



— ^ cosh 
2a 



\-)- 
a 



(7) 



Notice that with this parametrization the condition 
/ phdx = iVb is always satisfied so that a can be used 
as a fit parameter. In the left panel of Fig. [3 we have 
compared the bosonic density obtained using the self- 
consistent method with the Poschl- Teller approximation 
in Eq. |7J, from which we see that the approximation is 
indeed quite good. The right panel of the figure shows 
the gap soliton with four trapped fermionic levels in- 
side. Comparing these wavefunctions with the ones in 
Fig. IHI we see that, although different, they have the 
same symmetry properties with respect to the optical 
lattice. The fact that the bosonic density can be well ap- 
proximated by a Poschl- Teller potential allows to solve 
exactly the Schrodinger equation for the effective poten- 
tial Vef f = — 1x6/ 1 in terms of hypergeometric functions 
. The vanishing of bound state wavefunctions at large 
distances requires that the hypergeometric series must 
be reduced to a polynomial, this leading to the following 
analytical expression for the fermionic levels inside the 
condensate 



En = — 



^ (2n + l- ,Jl + 2a\xbf\Nb^ , n = 0, 1, n/. 

(8) 

The number of discrete levels nf in the gap soliton is then 
obtained as the largest integer n satisfying the inequality 



1) 



(9) 



In Fig. |S1 we compare the energy levels obtained with 
the self-consistent method with those obtained from Eq. 
©, from which we see that the agreement is quite good. 
In the right panel a similar comparison is made for the 
number of fermions Uf versus Xbf, this further con- 
firming the validity of our approximation. Notice that 
for Nblxbfla 3> 1 the energy levels © coincide with 
those obtained with the WKB approximation valid for 
n very large (in this limit Uf can be approximated as 
Uf ~ Nb\xbf\0'/'2)- From this we conclude that the 



number of fermionic levels at the bottom of the Fermi 
sea which can be trapped in a gap soliton can be con- 
trolled by changing the boson-fermion scattering length 
using Feshbach resonances. 




FIG. 8: Left panel. Energy levels of the fermionic bound 
states trapped inside an attractive condensate as a function 
of Xb/ • The continuous curves refer to numerical results ob- 
tained with the self-consistent method while the dotted lines 
to the analytical expression in Eq. ||SJ| with a = 0.66 and 
n — 0,1,..., 6, from right to the left, respectively. The top 
horizontal line represents the bottom of the fermionic band. 
Parameters are fixed as in Fig. |7| Right panel. Number of 
fermions entering the gap soliton as a function of Xbf- The 
continuous curve refers to Eq. @ with fitting parameter 
a = 0.66. To account for the existence of one bound state 
at small Xbf the curve has been shifted upward by a. All 
plotted quantities are in normalized units. 



CONCLUSIONS 

In conclusion, we have shown that localized states of 
BFM with attractive (repulsive) Bose-Fermi interactions 
can be viewed as a matter wave realization of quantum 
dots (antidots). The case of BFM in optical lattices has 
been investigated in detail and the existence of gap soli- 
tons has been shown. In particular, we showed that gap- 
solitons can trap a number of fermionic bound state lev- 
els inside both for repulsive and attractive boson-boson 
interactions. These solutions may have different symme- 
tries with respect to the OL and exist for a wide range 
of parameters. Gap solitons of BFM with OS and lA 
type are very stable under time evolution. We have 
shown that, as the boson-fermion interaction is changed 
the trapped fermionic levels inside the condensate may 
undergo level crossings or avoid crossings depending on 
the symmetry of the corresponding wavefunctions with 
respect to the optical lattice. This behavior is in agree- 
ment with a general theorem of quantum mechanics and 
shows that trapped fermions in BEC are in true quan- 
tum regime. The dependence of the number of bound 
states and of their energies on Nb and Xb / has been cal- 
culated both numerically and analytically. In particular, 
for attractive boson-boson interactions we approximated 
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the bosonic density with a Poschl- Teller potential and 
showed that a gap soliton in a shallow OL gives rise 
to a matter q-dot with filling number rif (x y^\xbf\Ni,. 
The possibility to change the boson-fermion interaction 
in a large interval of values (both negative to positive) by 
means of a Feshbach resonance, makes the filling depen- 
dence of a matter q-dot on Xbf experimentally feasible 
to check. To this regard we expect that as the (neg- 
ative) boson-fermion scattering length is decreased, the 
progressive filling of the gap soliton gives rise to a deple- 
tion of fermions in the OL which follow the above law 
and which could be in principle monitored by imaging 
techniques. To this regard we remark that the existence 
of gap-solitons in one dimensional Bose-Einstein conden- 
sates in optical lattices has been experimentally demon- 
strated in Ref . [lo| . We believe that similar experiments 
can be performed also for Bose-Fermi mixtures in OL 
and we expect that gap solitons will be found also in this 
case. 

Finally, we remark that the possibility to use gap soli- 
tons of BFM as q-bits protected from thermal dechoer- 
ence is very appealing and deserves further investigations. 
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